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Abstract 

We study the dynamics of a Hamiltonian system of N classical spins with infinite- 
range interaction. We present numerical results which confirm the existence of 
metaequilibrium Quasi Stationary States (QSS), characterized by non-Gaussian ve- 
locity distributions, anomalous diffusion, Levy walks and dynamical correlation in 
phase-space. We show that the Thermodynamic Limit (TL) and the Infinite-Time 
Limit (ITL) do not commute. Moreover, if the TL is taken before the ITL the sys- 
tem does not relax to the Boltzmann-Gibbs equilibrium, but remains in this new 
equilibrium state where nonextensive thermodynamics seems to apply. 
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1 Introduction 



Statistical thermodynamics is usually intended as the study of N-body sys- 
tems at equihbrium. However only a few textbooks [1] state clearly that the 
validity of equilibrium ensembles as models of thermodynamics is not auto- 
matically granted, but depends crucially on the nature of the Hamiltonian of 
the N-body system. In particular the very same basic postulate of equilibrium 
statistical mechanics, the famous Boltzmann principle S = k log W of micro- 
canonical ensemble, assumes that dynamics can be automatically (and kind of 
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easily) taken into account. However this is not always justified [2,3]. On the 
other hand, the Boltzmann-Gibbs canonical ensemble, is valid only for suffi- 
ciently short-range interactions and does not necessarily apply for example to 
gravitational or unscreened Coulombian fields for which the usually assumed 
entropy additivity postulate is not valid [4,5]. In general, a series of thermody- 
namic anomalies [6-12], which seem to escape a common general framework of 
understanding, has been observed. A few years ago, a generalized thermody- 
namics formalism based on a nonextensive entropy formula was proposed [13]. 
The latter has been encountering a large number of successful applications 
in far-from-equilibrium situations, as for example, to cite only a few cases 
among the most recent ones, in plasma physics [14], heavy-ion collisions [15], 
turbulence [16], discrete maps [17] and even in interdisciplinary fields such as 
bio-physics [18] and linguistics [19]. Such a formalism is the best candidate to 
be the general framework for a thermodynamics when long-range correlations 
or fractal structures in phase space are important and time evolution is not 
trivial, in other words when the dynamics plays a non trivial role [2,3]. In this 
paper we study the dynamics of relaxation to equilibrium in a Hamiltonian 
system of classical spins with infinite range interactions [20-23]. We show that, 
for some values of the initial energy and a class of off-equilibrium initial con- 
ditions, the systems does not relax to the Boltzmann-Gibbs equilibrium, but 
exhibits different equilibrium properties characterized by non-Gaussian veloc- 
ity distributions which can be fitted by the probability distribution functions 
(pdfs) of nonextensive thermodyamics[13]. The present study, together with 
the results presented in a paper now in press [24], provide the first indication 
that the generalized nonextensive thermodynamics can be a good candidate to 
explain some of the anomalies found in hamiltonian systems with long-range 
interactions. 



2 Dynamics and Thermodynamics of the HMF model 

The model, usually called Hamiltonian Mean Field (HMF), consists of N pla- 
nar classical spins interacting through an infinite-range potential [20]. The 
Hamiltonian is: 



where 6i is the ith angle and Pi the conjugate variable representing the rota- 
tional velocity. Note that the summation in V is extended to all couples of 
spins and not restricted to first neighbors. Following tradition, the coupling 
constant in the potential is divided by N. This makes H only formally exten- 
sive, i.e. V (X N when N — > oo[13,25], since the energy remains non-additive. 
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Fig. 1. The figure shows the initial time evolution for the Water Bag (WB) ini- 
tial state to the QSS for different system sizes. The straight full line indicates the 
canonical equilibrium temperature. 

that is the system cannot be trivially divided in two independent sub-systems. 
The model has an order parameter which is the magnetization M, i.e. the 
modulus of M = -i: I^ili where nij = [cos{9i), sin{9i)]. The canonical an- 
alytical solution of the model predicts a second-order phase transition from a 
low-energy ferromagnetic phase with magnetization M ~ 1, to a high-energy 
one, where the spins are homogeneously oriented on the unit circle and M ~ 0. 
The dependence of the energy density U = E/N on the temperature T, usually 
called the caloric curve, is given by [20,21] 




The critical point is at energy density Uc = 0.75 corresponding to a criti- 
cal temperature = 0.5 [20]. The dynamics of HMF can be investigated 
by starting the system with out-of-equilibrium initial conditions and integrat- 
ing numerically the equations of motion [21]. In particular in refs. [21,24] 
we have adopted Water Bag (WB) initial conditions, i.e. 9i = for all i 
[M = 1), and velocities uniformly distributed. In a special region of energy 
values (0.5 < U < Uc) the results of the simulations show, for a transient 
regime which depends on the system size, a disagreement with the canon- 
ical ensemble. In this region the dynamics is characterized by Levy walks 
and anomalous diffusion, while in correspondence the system shows a nega- 
tive specific heat [23] . Ensemble inequivalence and negative specific heat have 
also been found in self-gravitating systems [6,7], nuclei and atomic clusters 
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Fig. 2. For N=500 we show the time evolution of the temperature plateau corre- 
spondent to the Quasi Stationary State (QSS) calculated with different accuracies, 
reported in the plot. The initial time interval to = 100, shown in Fig.l, has been 
subtracted for clearness. 

[9-11], though in the present model such anomalies emerge as dynamical fea- 
tures [26,27]. In this paper we focus on a particular energy value belonging to 
the anomalous region, namely U — 0.69, and we study the time evolution of 
temperature, magnetization, velocity distributions. 

In Fig.l wc report the time evolution of 2 < A' > /N, a quantity that coincides 
with the temperature (< • > denotes time averages). The system, started with 
WB initial conditions, rapidly reaches a metastable Quasi-Stationary State 
(QSS) which does not coincide with the canonical prediction. In fact, after a 
short transient time, 2 < K > /N assumes a fixed value (the plateau in figure) 
corresponding to a N-dependcnt temperature Tqss{N) lower than the canon- 
ical prediction, also reported. In correspondence of this plateau one gets for 
the magnetization a value Mqss ~ 0. If we want to observe relaxation to the 
canonical equilibrium state with temperature Tcan = 0.476 and magnetization 
^can = 0.307, we have to wait for a time longer than that shown in fig.l, as 
shown for example in fig. 2 for the case N = 500. In ref. [24] the following 
scaling relations have been found: 

(i) the duration of the plateau, the lifetime of the QSS r, increases as r oc N; 



(u) TQssiN) 



0.380, a value obtained analytically as the metastable 



prolongation, at energies below Uc, of the high-energy solution (M = 0) as 
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Fig. 3. Time evolution of numerical velocity pdfs for N=1000 starting from WB 
and DWB initial conditions, panel a). For comparison we show also the equilibrium 
Gaussian curve (full line) in the QSS regime, panel b-c) (time t=1200,8000), and in 
the equilibrium one, panel d) (time t=500000). 



[Tqss{N)-T^] oc N-^/\ This implies that Mqss oc N-^/^ (see fig.l(c) of ref. 
[24]). 

These numerical results clearly indicate that 

• the two limits, the Infinite- Time Limit (ITL) t 00 and the Thermody- 
namic Limit (TL) N 00 do not commute] 

• if the TL is performed before the ITL, the system does not relax to the BG 
equilibrium and lives forever in the QSS. 

The robustness of the above results was checked in two different ways: 1) by 
adopting different initial conditions, as for example double water bag (DWB) 
initial conditions, and checking that we get the same QSS (see ref. [24] for 
details); 2) by changing the level of accuracy of the numerical integration, 
as shown in fig.2. We expect these results to be ubiquitous in nonextensive 
systems as conjectured in ref. [13]. 

In Fig. 3 we study the velocity pdfs. The initial velocity pdfs (WB or DWB 
initial conditions) quickly acquire and maintain during the entire duration 
of the metastable state a non-Gaussian shape. In Figs. 3(b) and (c) we see 
that the pdfs of the QSS do not change up to a time t = 8000 for a sys- 
tem with N = 1000 . The velocity pdf of the QSS is wider than a Gaussian 
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Fig. 4. Nonextensive theoretical pdf, with cut-off, (dashed curve) in comparison with 
the numerical ones (points) and the gaussian equilibrium prediction (full curve). In 
the inset we report the scaling of A = Pth — Pnum with N for those points indicated 
by the arrows. See text for further details. 

for small velocities, but shows a faster decrease for p > 1.2. The enhance- 
ment for velocities around p ~ 1 is consistent with the anomalous diffusion 
and the Levy walks observed in the QSS regime [22]. The following rapid 
decrease for p > 1.2 is due to conservation of total energy. From a dynam- 
ical point of view, the stability of the QSS velocity pdf can be explained 
by the fact that, for ^ oo, Mqss and thus the force on the spins 
Fi = {—M^sinOi + MyCosOi) — > 0. On the other hand, when N is finite, we 
have always a small random force,whose strength depends on N, which makes 
the system eventually evolve into the usual Maxwell-Boltzmann distribution 
after some time. We show this for N=1000 at time t=500000 in fig. 3(d). When 
this happens. Levy walks disappear and anomalous diffusion leaves place to 
Brownian diffusion [22]. 

We fit the non-Gaussian pdf in fig. 3(b) by using the one-particle prescription 
of the generalized thermodynamics (see ref. [24] for more details): 



pijp) 



1/(1-9) 



(3) 



This formula recovers the Maxwell-Boltzmann distribution for g = 1 and has 
been recently used to describe successfully turbulence [16] and non-Gaussian 
pdfs related to anomalous diffusion of Hydra cells in cellular aggregates [18]. 



6 



- 

- 

- 


U=0.69 QSS regime 












o ^-.^ 

/ ^ - ^ 






/ 












III 1 


1 



10^ 10^ 10'* 

N 



Fig. 5. Open points indicate the Largest Lyapunov Exponent (LLE) vs N for 
the QSS regime. An average over 10 runs is considered. Straight line indicates the 
theoretical prediction for LLE scaling, see text. 

In our case, the best fit is obtained by a curve with q = 7, T = 0.38 as 
shown in Fig. 4. The agreement between numerical results and theoretical 
curve improves with the size of the system. A finite-size scaling confirming the 
validity of the fit is reported in the inset, where A = Pth — Pnum, the difference 
between the theoretical points and the numerical ones, is shown to go to zero 
as a power of N (for four values of p). Since g > 3, the theoretical curve does 
not have a finite integral and therefore it needs to be truncated with a sharp 
cut-off (herein assumed to be discontinuous for simplicity) to make the total 
probability equal to one. It is however clear that, the fitting value q — 7 is 
only an effective nonextensive entropic index. Although similar non-Gaussian 
pdfs have been found previously, it was for dissipative systems [16], while this 
is the first evidence in a Hamiltonian system. 

In order to investigate deeper the dynamics of the plateaux observed in figs.l 

and 2, we have studied the Lyapunov exponent in the QSS regime. In fig. 5 
we show that, as expected, the Largest Lyapunov Exponent (LLE) tends to 
zero when N increases. The scaling behaviour of the LLE can be understood 
following the same argument already applied in ref.[21] in the overcritical re- 
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gion. It is known in fact that, when the Lyapunov can be estimated by the 
product of random matrices, the LLE scales with the power 2/3 of the pertur- 
bation [28,21]. In our case the perturbation is given by the magnetization, for 
which we have the scahng law oc N~^/^ in the QSS regime [24], thus we 
get the scaling LLE oc M^/^ oc (Ar-V3)i/3 ^ jy-^^. The latter is in perfect 
agreement with the numerical results as shown in the figure. Since LLE tends 
to zero as the system is increasingly large, one can safely say that mixing is 
increasingly slower and the observed anomalies in the relaxation process are 
naturally expected in the sense predicted by Krilov [29]. We note finally that, 
emergence of dynamical correlations and filamentary sticky structures in the 
//-space have also been observed in the QSS regime [24]. 



3 Conclusions 

In this paper we have studied a simple Hamiltonian system with long-range 
interaction. The dynamics of relaxation process is extremely rich and the sys- 
tem shows the existence of QSS different from the canonical equilibrium. These 
states satisfy the usual attributes of thermal equilibrium though they system- 
atically differ from what BG statistical mechanics has made famihar to us. 
Our results provide a first verification of nonextensive statistical mechanics 
[13] in long-range Hamiltonian systems, and illustrate the correctness of the 
criticism that Einstein developed in his celebrated 1910 paper [2] about the 
possible nonuniversality of Boltzmann thermostatistics, and the need of pro- 
viding a mechanical basis to it (including naturally time in the discussion) . 

Wc thank S. Ruffo for the suggestion of the LLE scaling in the QSS regime, 
and E.G.D. Cohen for stimulating discussions. 
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